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Abstract. - We consider liquid- vapor systems in finite volume V C K d at parameter values corresponding 
to phase coexistence and study droplet formation due to a fixed excess SN of particles above the ambient 
gas density. We identify a dimensionless parameter A ~ (5N) ( - d+1 ^ d /V and a universal value A c = 
A c (d), and show that a droplet of the dense phase occurs whenever A > A c , while, for A < A c , the excess 
is entirely absorbed into the gaseous background. When the droplet first forms, it comprises a non-trivial, 
universal fraction of excess particles. Similar reasoning applies to generic two-phase systems at phase 
coexistence including solid/gas — where the "droplet" is crystalline — and polymorphic systems. A sketch 
of a rigorous proof for the 2D Ising lattice gas is presented; generalizations are discussed heuristically. 



Introduction. - The thermodynamics of droplets in systems with phase coexistence has been well 
understood since the pioneering works [El-Q]. Recently, justifications of the classic results based on 
the first principles of statistical mechanics have been attempted — in both two and higher [fTo| — [l~3|] 
dimensions — and various thermodynamical predictions concerning macroscopic shapes have been rig- 
orously established. However, the formation and dissolution of equilibrium droplets is among the less 
well-studied areas in statistical mechanics. Indeed, most of the aforementioned analysis has focused 
on the situation implicitly assumed in the classical derivations; namely, that the scale of the droplet 
is comparable with the scale of the system. As is known [|7|,p|[l4]-|l6[|, this will not be the case when 
the parameter values are such that a droplet first forms. In this Letter, we underscore the region of 
the system parameters that is critical for the formation/dissolution of droplets. In particular, we iso- 
late the mechanism by which the low-density phase copes with an excess of particles and pinpoint 
the critical amount of extra particles needed to cause a droplet to appear. Surprisingly, at the point of 
droplet formation, only a certain fraction of the excess goes into the droplet; the rest is absorbed by 
the bulk. Moreover, apart from a natural rescaling to dimensionless parameters, all of the above can be 
described in terms of universal quantities independent of the system particulars and the temperature. 

In the last few years, there has been some interest in questions related to droplet formation and 
dissolution with purported applicability in diverse areas such as nuclear fragmentation JTt|— [T9|] and the 
stability of adatom islands on crystal surfaces [ |l"5| , |2(i| ] . Another issue, which is of practical significance 
in statistical mechanics, concerns the detection of first-order phase transitions by the study of small 
systems with fixed order parameter (magnetization) or fixed energy. Under these conditions, non- 
convexities appear in the finite-volume thermodynamic functions (which, of course, must vanish in 
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the thermodynamic limit), see j2l|— p4]] . Naturally, this suggests the formation of a droplet in a system 
with coexisting phases. Several studies have directly addressed the issues surrounding the appearance 
of droplets with intriguing reports on finite-size characteristics [|l5 , 20 , 2T,S,E^]. We believe that the 
results of this Letter may shed some light in these situations. 

Droplets in systems at phase coexistence. - We will start with some general considerations 
which bolster the claims of the first paragraph and, at the end of this Letter, describe the principal 
steps of a rigorous proof for the 2D Ising lattice gas. Although the natural setting for these problems 
is the canonical distribution in finite volume, intuition is often better developed in the context of finite 
subsystems using the language of the grandcanonical ensemble. Here the occurrence of droplets may 
be regarded as a problem in large deviation theory. This perspective will guide our heuristic analysis 
as it did in the proof for the 2D Ising system. 

Consider a generic liquid-vapor system. (In this Letter, we adopt, for concreteness, the language of 
the liquid-vapor transition. However, all considerations apply equally well to the formation/dissolution 
of an equlibrium crystal against a liquid or gaseous background.) First, suppose that the system is in the 
gaseous phase. According to a fluctuation-dissipation analysis, the local fluctuations for subsystems of 
volume V are then of the order s/xV, where, modulo constants, x is the isothermal compressibility. 
More precisely, the probability of observing a particle excess SN is given by 

, (SN) 
exp ' 



l 2xV J 



Now, when SN = 0(V xV), the above is just the leading-order asymptotic of a full-fledged Gaussian 
(central limit) distribution, which comes equipped with power-law corrections, etc. Moreover — in the 
single phase regime — the above leading order remains valid even for (SN) 2 ^ xV, provided that 
\SN\ <C PgV, where pc is the gas density. 

In the two-phase regime, small excesses can be again absorbed into background fluctuations but, 
in addition, a second mechanism exists through which the system can handle an excess of particles; 
namely, the formation of liquid-phase droplets. The minimal cost of a droplet of volume SV goes as 

expf-TwOHO^}, (2) 

where r\y denotes the (surface) free energy of an ideal-shape droplet of unit volume. For an isotropic 
system, t w = rSd, where r is the surface tension and Sd = 27r d / 2 /T(|) is the surface area of the 
unit sphere in R d . In general, t w is obtained by minimizing the Wulff functional [j3j]. As noted already 
in [Q-g|], by isoperimetric inequalities, scenarios involving multiple macroscopic droplets are far less 
likely. 

Now the number of excess particles 8N in a droplet of volume 8V is just 8N — (p^ — pc)5V, 
where pc and p^ are the ambient gas and liquid densities, respectively. Thus, comparing Eqs. (jj]) 
and (||), the droplet mechanism dominates when SN ^> QV d ^ d+1 \ where Q d+1 = (xTw) d (pL — 
pc) 1-11 , while the fluctuation mechanism dominates when SN <C QV d ^ d+1 \ We note that, from the 
perspective of rigorous analysis, significant progress has been made in the single-phase regime and in 
the above-mentioned extreme cases of the two-phase regime. For the Ising model, this was done for 
low temperatures in the exhaustive paper [Q], while [^|] extended this result throughout the coexistence 
region in the case d = 2. 

Previously, the crossover region SN ~ QV d '( d+1 ' has not received adequate attention. To study 
this region, we introduce the appropriate dimensionless parameter 
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and investigate finite (but large) size systems as A varies. The key to the whole picture is that, for the 
entire range of A, there is a forbidden interval of droplet sizes. To be precise, let us categorize droplets 
according to their surface area: We will say a droplet is of intermediate size if its surface area is large 
compared with log V but small compared with V^ 1 ^^ 1 ^. Droplets with surface areas outside this 
range will be called large and small as appropriate. We will show that, with overwhelming probability, 
there are no intermediate droplets. 

We begin with some observations: Suppose we specify the amount of excess which goes into 
intermediate and large scale droplets and fix the location of these droplets. Then, throughout the rest 
of the system, the fluctuations-dissipation result in Eq. ([[]) is valid, at least to leading order. Indeed, 
the only obstructions to Gaussian-type bulk fluctuations are: (1) The appearance of droplets beyond 
the logarithmic scale, which we have already separated for a special treatment, and (2) An exorbitant 
surface to volume ratio. Since SN « <~>y d /( d + 1 ) (i.e., A < oo), the second possibility does not 
occur. Let 

SN = 5N L + SNi + 5N S , (4) 

where SNl is the amount of excess particles in large droplets, and similarly for SN\ and SNs. By 
Eq. (|J), the distribution of the excess given in Eq. (Q) has a cost exp{— (SNs) 2 /(2kV)} for the fluctua- 
tion part. The cost of intermediate droplets will be of the order of their combined surface. If there are n 
such droplets, then isoperimetric reasoning forces us to pay at least exp{— TwC((57Vi)( £i ~ 1 )/ £i 7j 1/ ' d }, 
where C is a constant of order unity. On the other hand, if all of 8N\ were to go into small scale fluctu- 
ations, we would simply have to pay exp{— (SNs + SNi) 2 / (2xV)}. Comparing the two mechanisms 
we find, using SN S < SV d ^ d+1 '> and SNi < nQV d ^ d+1 \ that 

whenever n > 1. Hence, the probability of even a single droplet of the intermediate scale is utterly 
negligible. 

Having established the absence of intermediate-scale droplets, isoperimetric inequalities rule out 
the possibility of more than one large droplet. Thus, we are down to the simplest possible scenario: 
There is (at most) a single large droplet in the system — the cost of which is governed by Eq. (Q) — 
absorbing some of the excess, while the rest goes into background fluctuations — which are described 
by Eq. ([[]). Thus, the probability that the droplet contains the fraction A of the excess particles is, in 
the leading order, given by 



cxpi-r w (-^^) d *a(A)|, (6) 



where <E>a(A) is defined by 



$ A (A) = A" + A(l - A) 2 . (7) 



In particular, with overwhelming probability, the fraction of excess particles taken by the droplet 
corresponds to a value of A that minimizes $a(A). 

The result of a straightforward computation is that there is a constant A c , given by the expression 

i j , i d+l 

which separates two types of behavior: For A < A c , the unique global minimizer of $a(A) is A = 0, 
while for A > A c , the unique global minimum of $a occurs at a non-trivial value Aa > 0, see Fig. |jj 
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Fig. 1 - The graph of the universal function $a in d = 2. Here the parameter A represents the trial fraction of the 
excess that goes into the droplet; $a has the interpretation of afree energy function. In (a), A = 0.8 < A c , and 
the function is minimized by A = 0. In (b), A = 0.96 > A c and the function is minimized by a A = Aa > 2/3. 
The maximum that interdicts between A = and A = Aa presumably plays the role of a free energy barrier for 
the formation of the droplet. 



Moreover, the quantity Aa increases monotonically with A and the value of Aa at A = A c , denoted 
by A c , can be computed exactly; 

A < = dh- (9) 

In particular, we have Aa > A c for all A > A c . See Fig. ||. 

Let us interpret the results in the context of the canonical distribution: The region A < A c min- 
imized by A = is the remnant of the "phase" 5N <C QV d ^ d+1 ^\ the entire excess is taken up 
by background fluctuations. For A > A c , a large droplet occurs which absorbs the fraction Aa of 
the particle excess. Although this is obviously a precursor to the droplet-dominated "phase" (where 
SN ^ Qy d /i d + 1 ))^ the physics is somewhat different since a finite fraction of the excess — namely 
(1 — \<\)SN particles — is still handled by the background. We emphasize that Aa and A are related 
via a simple algebraic equation. The system-specific details and dependence on external parameters 
are encoded into the factor (p^ — pc) ( - 1 ~ d ^ d xTw from Eq. (||); the dimensionless parameter Aa is a 
universal function of A. 

We remark that in [14-16[], similar conclusions had been reached by various circuitous routes 



under the mantel of specialized assumptions or approximations. In this note, the exact formula has 
been derived on the basis of simple-minded droplet/fluctuation-dissipation arguments, all of which 
can be rigorously proved in at least one case. 

Mathematical results for 2D Ising model. - In the context of the two dimensional Ising lattice 
gas, the above reasoning has been elevated to the status of a mathematical theorem, which for conve- 
nience we state in the language of the equivalent spin system. Consider the square-lattice Ising model 
with the (formal) Hamiltonian 

w = - ax(J y> (10) 

where a x = ±1 and (x, y) denotes a nearest-neighbor pair. For each inverse temperature (3, let m* = 
m*(/3) be the spontaneous magnetization, \ — x(P) tne magnetic susceptibility, and tw = tw(/3) be 
the minimal value of the Wulff functional for droplets of unit volume. As is well known, m*(/3) > 0, 
< x(P) < oo and r w (/3) > once (3 > f3 = \ log(l + V2). 
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Fig. 2 - The graph of Aa, the fraction of the excess that goes into the droplet, as a function of A, the dimensionless 
rescaled parameter that measures the total excess. Here d = 2. Notice that Aa = for A < A c ~ 0.918, but as 
A I A c , Aa tends to A c = 2/3. The behavior of the system at A = A c has not been fully elucidated. 

Consider now an L x L square in Z 2 denoted by and let Ml = J2 x eA L a % ^ e tne ove rall 
magnetization in A^. Let vl > be such that m*| Ax,| — 2m*VL is an allowed value of Ml for all L. 
Let Pl'v l t> e tne canonical distribution on A^ with plus boundary conditions, inverse temperature (3, 
and Ml fixed to the value m*|Ai| — 2m* Vl- In the present setting, the parameter A in Eq. (j^) becomes 

I *\2 3/2 

(m ) v T ' 
A = 2i-^- lim t-^j, (11) 

xm l^oo \Al\ 

where we presume that the limit exists. 

In Ising systems, a convenient description for the spin configurations is in terms of their Peierls' 
contours, i.e., the lines separating spins of opposite type; see Fig. [|. Moreover, in the present context, 
the boundaries of droplets are exactly these contour lines. Our first claim concerns the absence of 
contours of intermediate size, regardless of the value of A. 

Theorem I. Let (3 > f3 c and suppose that the limit in Eq. f |77| ) exists with A G (0, oo). Let Al be the 
event that there is no contour T in Al with 

KlogL < diamT < -^£ 2/3 . (12) 
K 

If K — K(f3) is sufficiently large, then 

lim P+ : l(A L ) = l. (13) 

We remark that this is the rigorous (albeit 2D-Ising specific) analogue of our general argument 
in Eq. ([}])-Eq. (^). The above theorem is far and away the most difficult part of the mathematical 
analysis. Notwithstanding, the flow of the proof parallels closely the derivation that was given here. 
The reader may have noticed that, in the present derivation, we have not used in any obvious way the 
lower bound defining the scale of intermediate droplets. The issue is somewhat delicate since some 
log-scale droplets will naturally emerge from the background fluctuations. The constant K must be 
chosen large in order to enforce the distinction between "natural" and "unnatural" as well as to control 
the translation entropy of the purported intermediate droplets. 

Our next goal is to specify the typical configurations in measure Pl'v l depending on the value A 
as compared with A c . Let Bl, k be the event that there is no contour T in Al with diam T > K log L. 
Furthermore, let Cl,k be the event that there is one contour Tq with diamTo > j^L 2 / 3 while all other 
contours T in Al satisfy diam T < K log L. 
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Fig. 3 - Peierls' contours in a 2D-Ising spin configuration. In the lattice-gas language, we regard minus spins as 
particles and plus spins as vacancies. Peierls' contours are then the interfaces separating high and low density 
regions at the microscopic scale. 



Theorem II. Let (3 > (3 C and suppose that the limit in Eq. ( |7"/| ) exists with A G (0, oo). Let A c be as 
in Eq. (js]) and, for A > A c , let Aa be the unique minimizer o/$a(A). 

(1) If A < A c and K — K{j3) is sufficiently large, then 

lim P££(B L , K ) = 1. (14) 

(2) If A > A c and K = K{(3) is sufficiently large, then 

lim Pt;i{C L , K ) = l. (15) 

Moreover, with probability approaching one, the unique "large" contour Tq has volume (\a+o(1))vl 
and its shape asymptotically optimizes the surface-energy (Wulff) functional for the given volume. 

Theorems I and II completely classify the behavior of the Ising system for all A A c . We 
emphasize that the situation at A = A c has not been fully clarified. What can be ruled out, according 
to Theorem I, is the possibility of a complicated scenario involving intermediate size droplets on a 
multitude of scales. Indeed, when A = A c , in (almost) every configuration we must have either a 
single large droplet or no droplet at all; i.e., the outcome must mimic the case A > A c or A < A c . 
It is conceivable that one outcome dominates all configurations or that both outcomes are possible 
depending on auxiliary conditions. 

Our last statement concerns the decay of the probability (in the grandcanonical distribution) that the 
overall magnetization takes value m*|Ai| — 2m* v^. Let be the (Gibbs) probability distribution 
on spins in Al with plus boundary condition and inverse temperature (3. 

Theorem III. Let (3 > (3 C and suppose that the limit in Eq. ( |7i| ) exists with A € (0, oo). Introduce the 
shorthand = ^'^(AIi = m*|A^| — 2m*vi J ). Then 

lim v L \ogp L = -t w inf $a(A). (16) 
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Sketch of the proofs. - Here we outline the steps necessary to prove the above theorems. As 
already noted, first we reduce the problem to the study of large-deviation properties of the "grand- 
canonical" distribution J 1 ^ '* using the relation 

Pl$M) = ^(AMl = m*|A L | - 2m*v L ), (17) 
valid for all events A. The next technical step is then a proof of the large-deviation lower bound 

PL > exp{-r Wx /^I($A(A) + e L )}, (18) 

which is produced by forcing in a contour of the appropriate size and evaluating the contributions from 
surface tension and bulk fluctuations. Here — * as L — ► oo, uniformly in A. A comparison with 
this lower bound then shows that, with overwhelming probability in P^'„ , the total surface area of 
contours T with diam V > K log L is at most of order y/vZ while their combined volume is at most 
of order vl- This puts us in a position to carry out the argument in Eq. (|j), which ultimately leads 
to the proof of Theorem I. Having eliminated the intermediate contours, we are down to the scenario 
with at most one large contour. Optimizing over the contour volume/shape proves Theorem II and also 
produces a large-deviation upper bound, which completes the proof of Theorem III. 

Theorems I-III pretty much tell the story for this particular case. Complete proofs and additional 
details will appear elsewhere [^7[], 

* * * 
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